1. IntroductionManipulation of atomic Bose–Einstein condensate (BEC) is very important in many applications, such as atom laser,[1,2] matter-wave interferometry,[3,4] atomic holography,[5–7] and lithography.[8,9] In particular, to obtain a matter wave with specific de Broglie wavelength, accelerating or decelerating of BEC is regularly needed. Very basic mechanical knowledge states that a free wave packet can only keep still or move at a constant velocity, to accelerate it an external potential should be applied. {Along} this line of thought, many accelerating schemes have been proposed, for example, using of gradient magnetic field,[10] time-dependent harmonic traps,[11] frequency-chirped optical lattice,[12] spatial array of quadrupole traps,[12] and modulated optical fields.[13]
Moreover, it is found that certain kind of accelerating waves can also be realized without any external potential applied. In the year 1979, Berry and Balazs theoretically investigated the free space Schrödinger equation and showed that an Airy wave packet can undergo an accelerating motion in free space.[14] This miraculous property of the Airy wave packet originates from its non-integrability (in other words, an Airy wave packet carries infinite amount of matter and energy, therefore is not physically realizable). After Berry and Balazs’s pioneering work, extensive theoretical investigations were done on such kind of waves (for a review see Ref. [15]). And it was pointed out that an Airy wave packet can be made integrable (thus physical realizable) and at the same time holding its main feature by truncating to a finite range.[16] In 2007, such self-accelerating truncated Airy wave packets were experimentally demonstrated in an optical system.[17] The idea was soon introduced back to the field of matter wave. Accelerating electron beams have already experimentally achieved.[18] Recently, the generation of accelerating BEC Airy wave packets by using amplitude or phase imprinting techniques was also proposed.[19]
Now, the question is besides these two ways, are there any other methods to accelerate a BEC wave packet? For atomic BEC, the inter-atom interaction can be tuned by the Feshbach resonance technique,[20,21] and is also an important and widely used tool for the manipulation of BEC.[22–31] But it is generally believed to have no effect on accelerating of a matter wave packet,[11] since the interaction between two particles usually conserves the total momentum. In this article, at first, we theoretically show that this belief is valid only for spatially homogeneous interactions (regardless of the particular interaction form and wave packet shape), while in the case of spatially inhomogeneous interactions it is broken. As an example, dynamic of BEC under spatially linear modulated contact interaction is studied. We show that an initial steady BEC wave packet will undergo an accelerating motion under such an interaction. But the acceleration is a time varying one. To overcome this blemish, we also propose an engineering scheme of the interaction to achieve a constantly accelerating matter wave packet. Numerical results demonstrate that this scheme works well, and at the same time this scheme also suppresses profile varying of a wave packet during its evolution, thus produces a profile-keeping accelerating matter wave packet.
The rest of this article is organized as follows. In Section 2, the effect of the inter-atom interaction on accelerating dynamic of a BEC wave packet is studied theoretically. In Section 3, accelerating motion of a BEC wave packet under spatially linear modulated contact interaction is studied numerically. In this part, a scheme to realize a constantly accelerating profile-keeping wave packet is also proposed. In Section 4, the main results of this work are summarized.
2. TheoryTo simplify the discussion, here we consider a system of BEC which is tightly trapped along the y and z directions with an external potential (
), yet can freely move along the x direction. Such a system can be reduced to one dimension,[32,33] and under the mean field theory the dynamic can be well described by the following Gross–Petaevskii equation:
where
ħ is the reduced Planck constant,
px = −i
ħ∂/
∂ x is the momentum operator along
x direction,
m is the mass of the condensated atom,
N is the total number of atoms contained in the condensate,
V(
x,
x′) is the interaction of two atoms at positions
x and
x′ respectively, and xde
ψ is a wave function describing the macroscopic quantum state of the condensate which is normalized to 1, i.e.,
.
By using Gross–Petaevskii equation (1), the equations governing the evolution of the mean values of the BEC wave packet position and momentum, ⟨x⟩ and ⟨p⟩, can be derived as (Ehrenfest’s theorem[34,35])
with
Veff(
x,
t) =
N∫
V(
x,
x′)|
ψ(
x′,
t)|
2d
x′ being an effective potential produced by the inter-atom interaction. Comparing with the classical correspondence, it is clear that
v and
a are the moving speed and acceleration of the BEC wave packet, respectively. If
a = 0, the wave packet will stay still or moving at a constant speed. However, if
a ≠ 0, the wave packet will move in an accelerating manner. Using commutation relation [
px,
f(
x)] = −i
ħd
f(
x)/d
x, equation (
3) can be further calculated as follows:
Since there are two space variables
x and
x′ (positions of two interacting atoms) in the expression of acceleration, it would be convenient to introduce new variables
In terms of
X and
xc, acceleration (
4) reads
which can be split into two parts
with
Because we are considering the interaction between two identical particles, switching their positions will make no difference, i.e., regardless of the special form of interaction we always have V(X,xc) = V(−X,xc). This is to say, the integrand in Eq. (9) is an odd function with respect to variable X, thus the value of a1 will always be zero. It will have no contribution to the acceleration of the BEC wave packet. When the interaction is spatially homogeneous, that is, the interaction potential V has no relation with {variable} xc, or in other words, ∂ V(X,xc)/∂ xc equals zero. Then acceleration a2 will also be zero. So, a spatially homogeneous interaction will have no accelerating effect on the BEC wave packet regardless of the specific form of the interaction (contact, dipole–dipole, or some others). This is the usually well known case.
Here it should be pointed out that for a self-accelerating Airy wave packet, because its wave function is not square integrable, consequently there does not exist proper definitions of ⟨x⟩ and ⟨p⟩, and the above derivation fails fundamentally. And for the square integrable truncated Airy wave packet, although the wave packet seems to undergo an accelerating motion, its mean position ⟨x⟩ will have no change, and its acceleration a will keep zero during its propagation (see Fig. 1). That is, the truncated Airy wave packet in fact does not accelerate, therefore we call the motion a seemingly accelerating motion.
However, when the interaction is spatially inhomogeneous, a2 can have a non-zero value, and the condensate will be accelerated by the interaction. One such example can be easily fabricated
where
δ(·) is the usual Dirac delta function (the interaction is a commonly considered contact interaction), and
g(
x) is a space-dependent interaction strength. Bu using the Feshbach resonance technique,
[20,21] such an interaction can be realized by applying a spatially modulated magnetic field.
[23–31] Here we study a simple linear form of the interaction modulation
where the parameter
g0 = 2
ħω⊥as/
l0 (
as is the s-wave scattering length of the condensate atom, and
l0 is the length scale of the interaction gradient) is an interaction strength scale, and
xg = 0 defines the point at which the interaction strength is tuned to 0.
Firstly, to roughly get some physical insight, we consider a Gaussian wave packet
with
x0 and
σ0 being the center position and width of the wave packet, respectively. For such a wave packet, the acceleration can be calculated analytically, and it is
which is obviously not zero, indicating that the BEC wave packet can be accelerated by the interaction induced effective potential. Considering a
39K BEC with 10
6 atoms
[36] transversely trapped by a harmonic trap with
ω⊥ = 1000 Hz and having length
σ0 = 5 μm along
x-direction, if the s-wave scattering length
as is tuned from 0 to −25
a0 within length
l0 = 1 cm (for
39K,
as = −33
a0(1−52/(
B − 402.4)) with
B in units of G,
[37] so this can be achieved by nonlinearly tuning the magnetic field from 454.4 G to 604.4 G within length of 1 cm, the corresponding magnetic gradient is in the order of magnitude of 150 G/cm, which is ordinarily used in BEC experiments
[38–40]), the acceleration is estimated to be about ∼ 0.08 m/s. This value can be greatly increased if a large atom number (for example,
N = 10
8[41,42]) BEC or a magnetic field with large gradient (for example, 200 T/m
[43]) is used.
It should be noted that unlike accelerated by a linear external potential, here the acceleration not only depends on slope g0 of the interaction potential, but is also related to width σ0 of the wave packet. Under the same interaction, a narrower (smaller value of σ0) wave packet will have a larger acceleration. Moreover, if the evolution of such an initial wave packet is considered, as the width of the wave packet will usually change during its evolution due to the dispersion and interaction, the interaction induced acceleration will also be a time varying one (this will be shown in the next section in detail numerically).
Although equation (14) comes from a Gaussian wave packet assumption, it is valid qualitatively for wave packets with other mathematical forms, such as sech which is the mathematical form of a bright soliton.[22,44,45] This can be explained as follows. Because of the normalizing restriction (i.e., conservation of particles number), a smaller value of wave packet width σ0 will produce a sharper effective potential, and it will induce a larger force and hence acceleration on the wave packet.
3. Numerical resultsTo study the dynamic of such a system, one substitutes the specific interaction potential (11) into Eq. (1), and immediately obtains the following nonlinear Schrödinger equation:
It is usually hard to solve this equation analytically, so here we solve it numerically using the extensively used operator splitting and fast Fourier transform technique.
[46,47] And in the calculations, natural units
ħ =
m = 1 will be used.
In Fig. 2, we take function (13) with parameters x0 = 0 and σ0 = 1 as the initial wave packets and show their evolutions under different forms of interaction. In Figs. 2(a1)–2(a4) and 2(b1)–2(b4), the interactions are set to Ng(x) = 0.3(x + 10) and Ng(x) = −0.3(x + 10), respectively. It is clear that in these two cases, the BEC wave packet undergoes an accelerating motion, the ⟨x⟩–t relation shows a curved shape and the acceleration a has a non-zero value. At the very beginning time, the shape of the wave packet has not changed much, thus the value of acceleration a agrees well with the value given by Eq. (14), a = − 0.042 and 0.042 for the two cases, respectively. In Figs. 2(a1)–2(a4), the interaction has a net repulsive effect in the condensate occupied space. Both the interaction and dispersion effects tend to broad the wave packet during its evolution, thus according to Eq. (14), the magnitude of acceleration decreases. While, in Figs. 2(b1)–2(b4), the interaction has a net attractive effect in the BEC occupied space. In such a case, the attractive interaction and the dispersion effects compete with each other, as a result the width of the wave packet experiences an oscillation during the evolution, and consequently the acceleration oscillates as well. To show this relation between the wave packet width and acceleration, in Fig. 2(b4) we also plot the full width at half maximum (FWHM) of the wave packet. The correlation between a and FWHM can be seen clearly in this figure. For comparison, in Figs. 2(c1)–2(c4) we plot the evolution of the same initial wave packet under a spatially homogeneous interaction Ng(x) = −3. We see that although the wave packet breathes during its evolution, it does not accelerate. Its mean position ⟨x⟩ and acceleration a keep 0 all the time.
To demonstrate that the acceleration in Figs. 2(a) and 2(b) is really caused by the interaction, we plot the corresponding kinetic and interaction energies separately in Fig. 3. From the figure, we see that the increasing of the kinetic energy is always accompanied by a decreasing of the interaction energy. This fact clearly indicates that the acceleration of the wave packet is due to the inter-atom interaction. The oscillation of Ek,i in Fig. 3(b) is also due to breathing of the wave packet.
From Figs. 2(a4) and 2(b4), it also can been seen that the accelerations are time varying during the evolution. That is, the motions of the wave packets are not usually interested simple constantly accelerating ones. This blemish can be overcome by engineering the interaction strength factor g0 over time using the Feshbach resonance technique.[20,21] To achieve a constantly accelerating motion with acceleration a = ac, according to Eq. (14), initially one should set the interaction strength to be
and at time
t one needs to tune the interaction strength to
with
In Fig. 4, we plot the engineered g0(t) and the corresponding dynamical evolution of a wave packet (13) under the attractive interaction case. From the figure, we see that under the engineered interaction the acceleration keeps its initial value, and the mean position increases in a perfect parabolic manner. Thus, a constantly accelerating wave packet is achieved as expected.
Moreover, comparing Fig. 4(b) with Fig. 2(b2), we can roughly notice that the interaction engineering can also suppress shape varying of the wave packet during its evolution. This can be seen more clearly from Fig. 5 where we explicitly compare the FWHM of the same wave packet under constant and engineered interactions during the evolution. From this figure, we see that after t = 2 the width of the wave packet changes very slightly under the engineered interaction. Thus, the achieved constantly accelerating wave packet also has a nearly shape-keeping feature. These manners make the system potentially useful in the field of creating atom laser[1,2] and atomic analogy of light bullet.[44,48]
Although for the repulsive case, one can also get a constantly accelerating BEC by engineering g0 in a similar way (see Fig. 6), the engineering will at the same time speed up the broadening of the wave packet (to keep the acceleration constant, the broadening of the wave packet needs a stronger interaction to balance it, and a stronger repulsive interaction will make the wave packet broad more quickly), see Fig. 7 where we compare the FWHM of wave packets during evolution under constant and engineered interactions, i.e., the wave packets in Figs. 2(a2) and 6(b). The broadening of the wave packet will greatly reduce the intensity of the matter wave, therefore may make it less significant in practice.
4. ConclusionWe discussed the role of interaction in accelerating a BEC wave packet. We showed that if the interaction between two atoms is spatially homogeneous (or in other words, only depends on their separating distance), it will have no acceleration effects on the atomic BEC. However, if this is not the case, the interaction may cause an accelerating motion of the condensate. One example of such interaction induced accelerating BEC was fabricated. We showed that if we linearly modulate the usual contact interaction in space, an initial steady BEC wave packet will undergo an accelerating motion with time varying acceleration. We also proposed that a constantly accelerating BEC wave packet can be obtained by timely engineering the interaction strength. At last, the analysis here can also be straightforwardly extended to the field of nonlinear fiber optics where a light pulse obeys the similar nonlinear Schrödinger equation.[49]